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Abstract. A microscopic theory of superfluid drag in a two-component Bose 
gas is developed. The drag factor is shown to be proportional to the square 
root of the gas parameter. Basing on the similarity between the "drag force" 
and vector potential of magnetic field we propose how the superfluid drag can 
be used to detect a Berry phase in a Bose counterpart of the Josephson charge 
qubit. We consider a system in which the drag component, situated inside the 
drive component, is confined in two half-ring traps separated by two Josephson 
barriers. Under cyclic adiabatic change of Josephson coupling parameters and an 
asymmetry of two traps a Berry phase is generated. This phase can be observed 
though measurements of relative number of atoms in two traps. The Berry phase 
depends on the drag factor and its detection can be used for determining the value 
of the drag. 
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1. Introduction 

The interest in studying of macroscopic quantum phenomena in Josephson 
superconductive devices has renewed recently in connection with the possibility of 
using such devices for performing quantum logic operations p^. Related problem - 
Josephson-type effects in Bose systems was received considerable attention [5] after 
experimental realization of Bosc-Einstein condensation in ultracold atomic gases. It 
was shown in SJ that Josephson Bose system can also be used for quantum computing. 
In [3] a simple model of a Bose counterpart of the Josephson charge qubit has been 
formulated. 

The supercurrent in superconductors is coupled to the vector potential of 
electromagnetic fields that allows to perform quantum logic operations by external 
magnetic fields. External magnetic fields couple with Bose atoms only through Zeeman 
and hyperfine interactions. Therefore, in difference with the case of superconductors, 
magnetic field cannot be used to control supercurrents in uncharged Bose systems. 

In this paper we consider the effect that can be used for quantum state engineering 
in Bose Josephson systems in the same manner as magnetic field in superconductive 
Josephson qubits. This effect is a nondissipative drag of supercurrent in a two- 
component superfluid system. 

Nondissipative drag is caused by interaction between the atoms and consists in 
a redistribution of a supercurrent between two superfluid components. For the first 
time the nondissipative drag was considered by Andreev and Bashkin [J. In 0] 
a three-velocity hydrodynamic model of a 3 He- 4 He mixture was developed. It was 
shown that superfluid behavior of such a system can be described under accounting 
the " drag" term in the free energy. This term is proportional to the scalar product of 
the superfluid velocities of two components times the difference between the effective 
and the bare masses of 3 He atoms. A microscopic theory of the nondissipative drag 
for a special case of a bilayer system of charged bosons was developed in Refs. |SJ|S]. 

Nondissitive drag is a fundamental property of systems with a macroscopic 
quantum coherence and takes place also in superconductors , quantum Hall systems 
[B] and quantum wires 9 . But for the system of uncharged bosons this effect is 
especially important since the "drag force" plays the role similar to the role of the 
vector potential of a magnetic field in superconductors 4 . 

An important problem of quantum computation is to minimize the influence of 
random noise coming from the environment. One of the ways of solutions of this 
problem consists in using of geometrical phases to implement quantum gates . As 
was shown in [IT| a geometrical Berry phase can be generated in a SQUID loop with 
two Josephson junctions in an external magnetic field. In this paper we propose a 
model of Bose qubit that may demonstrate analogous effect. 

In Sec. El the theory of nondissipative drag in two-component Bose gases is 
developed. In Sec. Olthe model of Bose qubit is formulated. The scheme of generating 
of the Berry phase is proposed in Sec. ^ Conclusions are given in Sec. El 

2. Nondissipative drag of superflow in a two-component Bose system 

In this section we develop a microscopic theory of the nondissipative drag in a two- 
component atomic Bose system with point two-particle interaction. 

Let us consider an infinite uniform two-component atomic Bose gas in a Bosc- 
Einstein condensed state. In general case the densities of atoms in each component 
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are different from one another (ni ^ n 2 ), the atoms of each components have different 
masses (mi 7^ 7712) and the interaction between atoms is described by three different 
scattering iengths (an 7^ 0,22 7^ 012)- The Hamiltonian of the system can be presented 
in the form 



i=l,2 



i,»'=l,2 



where 



Ei = I d 3 r— [W+(r)]W,(r) 



is the kinetic energy, 



d s r*+(r)*J(r) 7 «,* i ,(r)*i(r) 



(1) 



(2) 



(3) 



is the energy of interaction, 7,^ = 4ttH 2 au/rrii and 712 = 2TtTi 2 (m,i + 7712)012/(77117712) 
are the interaction parameters. 

For further analysis it is convenient to use the density and phase operator 
approach (see, for instance, The approach is based on the following 

representation for the Bose field operators 

exp [i<pi{r) + i(fii(r)] v /n l + fn(r), (4) 

(5) 



*+(r) = y/m + hi(r) exp [-i</?i(r) - i&(r)] , 



where n, and ^ are the density and phase fluctuation operators, <fii(r) is the c- 
number term of the phase operator, which are connected with the superfluid velocity 
by common relation Vi = hWfi/mi. In what follows we specify the case of superfluid 
velocities independent of r. 

Substituting Eqs. J3J, (JSJ into Eq. Q and expanding it in series in powers of hi 
and V</?i we present the Hamiltonian of the system in the following form 

H = H + H 2 + ... (6) 

In jnj the term 



He 



1 2 , 7it 2 

-mmjVj + — - ;Ltj?l, 



712771772 



(7) 



does not contain the operator part (V is the volume of the system). The chemical 
potentials are determined from the minimization condition of the Hamiltonian Hq: 

-rn,v? + 7iin, + 712^3-, - /i, = 0. (8) 

Fulfillment of Eq. ((HJ results in vanishing of linear in the density fluctuation operators 
terms in the Hamiltonian. The terms linear in the phase fluctuation operators are 
absent in the Hamiltonian due to continuity conditions. 

The quadratic in and hi operator term in H reads as 

2 



Hi 



/*(? 



2m,- 



(VrW(r))' 

477,, 



+ 77. t (v^(r))' 



+-| i (n;(r)V&(r) + [V&(r)]n;(r)) 
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^([Vni(r)]V&(r) - [Vft(r)]Vn<(r)) + y ("*(')) ^ [ + 7i2"i(r)n 2 (r) 



(9) 



The quadratic part of the Hamiltonian determines the spectrum of the elementary 
excitations. Hereafter we will neglect the higher order terms in the Hamiltonian JBJ. 
These terms describe the scattering of the quasiparticles and they can be omitted if 
the temperature is much smaller than the temperature of Bose-Einstein condensation. 

Let us rewrite the quadratic part of the Hamiltonian in terms of the operators of 
creation and annihilation of the elementary excitations. As the first step we use the 
substitution 



'•Pi 



(10) 

(11) 



where operators bf, bi satisfy the Bose commutative relations. Here en- — H 2 k 2 /2mi 
is the spectrum of free atoms, and Eik = y 1 eiki^ik + ^Jan-i) is the spectrum of the 
elementary excitations at 712 = and Vj = 0. The substitution reduces the 

Hamiltonian Q to the form, quadratic in operators bf , 6$: 



Ho 



E 



St(k) 6+(k)6,-(k) 



Here 



and 



+ E-? fc K( k ) & 2(k) + &i(k)6 2 (-k) + h. 

k 

£i(k) = E lk + Kkv, 



(12) 
(13) 

r r - (14) 

The Hamiltonian 112|l contains non-diagonal in Bose creation and annihilation 
operator terms and it can be diagonalized using the standard procedure of u-v 
transformation I 1 . The result is 



9k = 712 



eife£2feni^2 



Ho 



E 



J2 fA(k)(/3+(k)/3 A (k) + J] -\ e« 



1=1,2 



(15) 



where P\(k) + and f3\(k) are the operators of creation and annihilation of elementary 
excitations. 

The energies £a (k) are found as solutions of the equation 
A £1 B 



det 



B 



A + £I 



= 0, 



where 



/£i(k) 





9k 








^i(-k) 





9k 


9k 





f 2 (k) 





\ 


fjk 





&(-k) 



(16) 



(17) 
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(18) 



and I is the identity matrix. Having the spectra of elementary excitations one can 
obtain the relation between the supercurrents in each component are their superfiuid 
velocities. 

We will find the current densities using the relation 



. _ 1 8F 



(19) 



where F is the free energy of the system. Here the quantity is defined as density of 
the mass current. 

The free energy of the system, described by the Hamiltonian (JBJ) , is given by the 
formula 



F = H { 



o + jE 



J2 &(k) - £ e ik 

X=a,0 i=l,2 



1 — exp 



gA(k)\ 



(20) 



k \=a,/3 

The second term in (I20|l is the energy of the zero-point fluctuations and the third 
term is the standard temperature dependent part of the free energy for the gas of 
noninteracting elementary excitations. 

To proceed further we consider the case of small superfiuid velocities, much 
smaller than the velocities of elementary excitations, and seek for the free energy 
as series in v,;. 

At Vi = v-2 = the equation l|16fl is easily solved and the spectra are found to be 

\ 1/2 



E, 



a{f3)k 



E 2k ± 



+ 47j J 2 nin 2 eifce 2 fc 



(21) 



As required in the procedure ^1] we take positive valued solutions of Eq. (jTB|) . We 
should emphasize that the energies l|21|) should be real valued. This requirement yields 
the common condition of stability of a two-component system: jf 2 ^ 7n 722- If this 
condition were not fulfilled spatial separation of two components (at positive 712) or 
a collapse (at negative 712) would take place. 

At nonzero superfiuid velocities we present the solutions of Eq. (|16fl as series in 
Vj, neglecting the terms higher than quadratic: 



£ a (k) = E ak 



-nk Vl 1 



p2 

E lk ~ E 2k 



E lk 



— Tikv? 
2 



27^ 2 nin 2 ei fe e 2 fe 

(3^ fc + E} k ) 



EnU. I E. 



E lk, 



ft 2 (kv! — kv 2 ) 



(?2 p2 
E lk - E 2k 
p2 _ p2 

^ak ^/3fe / 



(22) 



^(k) = E i 



0k 



-fckvr 1 



Elk 



E 2k 



27f 2 nin 2 eifee 2 fc ( E~ 



El 
3E jk 



F 2 

^Pk 



-?Lkv 2 1 



E ak ( E ik- E j k y 



h 2 (kvi — W2) 



p2 rp2 
E lk - E 2k 



Ei 



E 2 



(23) 
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Note that at Vj = v 2 = v the spectra Q22[l. (|23[1 are reduced to common expressions 
for the energies of quasiparticles in a moving condensate: £ a /p\(k) = E a fp\ k + frkv. 

Using Eqs. H2U|) . il22[l an d we obtain the following expression for the free 
energy 



F = F 



V 
1 



(Pi 



where Fq does not depend on v, 
quantities 



Pnl)v\ + (P2 
In fify 0s = 



P«2)V2 - Pdr(vi - V 2 ) 2 



(24) 



Pnl = 


2mi 


k 


dN ak 


h dN pk 




~ S 2fc 


(dN ak 


dN pk 




dE ak 


dEp k 






K dE ak 


dE 0k 


Pn2 = 


2m 2 
~ 3V 


k 


dN ak 
dE ak 


dN pk 
dEp k 




~ E 2k 
~ E lk 


( dN ak 
V dE ak 


dNp k 
dEp k 



rrijUi are the mass densities, the 

(25) 
(26) 



describe the thermal reduction of the superfluid densities, and the quantity 



Pdr = — V TO 1™2 2_^ 
'* k 



2 i ^3/2 

7i2 n i"2 (eifee 2 fc) 



E nk E, 



k^Bk 



1 + N ak + N, 



0k 



N, 



0k 



{E ak + E pk f (E ak - E pk y 



2E ak E 0k fdN ak dN pk 



{Elk- E}k) 



\dE a k 



dE 



3k 



(27) 



which we call the "drag density," yields the value of redistribution of the superfluid 
densities between components. In Eqs. I|25 |) -l|2"7 )l N a rm k — [exp(E a rm k /T) — is 
the Bose distribution function. 

Using Eqs. i|19|) . jUt we arrive to the following expressions for the supercurrents 

jl = (Pi - Pnl - Pdr)^! + PdrV 2 , (28) 
(P2 - Pnl - Pdr)v 2 + PdrVl. (29) 



J 2 



One can see that at nonzero pd r the current of one component contains the term 
proportional to the superfluid velocity of the other component. It means that there 
is a transfer of motion between the components, or the drag effect. Since pdr oc jf 2 
the effect occurs only at nonzero interaction between the components and it does not 
depend on the sign of this interaction. 

Eq. (|27|l is the main result of this section. This equation yields the value of 
the drag for a general case of two-component Bose system with components having 
different densities, different masses of atoms, different interaction parameters, and for 
zero as well as for nonzero temperatures. 

To estimate the absolute value of the drag we, for simplicity, specify the case of 
mi = m 2 , that is realized when two components are two hyperfine states of the same 
atoms. 

At mi = m 2 = m and T = Eq. I|27|l is reduced to the form 

-yf 2 mn 2 iy(e)e 1/2 



4m . 
Pdr = — I de 



o i/(e + w 1 )(e + w 2 ) (V e + w\ + y/e + w 2 ) 



where 



wi{2) = 7n™i + 722«-2 ± y (7n"-i - 722"2) 2 + ^l 2 n]_n 2 



(30) 



(31) 
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and 

vie) = —= 5-Ve 

V2n 2 h 3 

is the density of states of elementary excitations. 

The integral in l|30|) is evaluated analytically. To present the answer in a compact 
form it is convenient to introduce dimensionless parameters 



,2 



rj = ^ 12 and n = . — — — . 

711722 V n 2722 

Using these notations we have 



Pdr = ^/pip2\fn 1 al 1 n 2 al 2 r] 2 F(K,T7), (32) 

F( K)?7 ) = -^=(« ; + K- 1 +3v / r^ 
45V27T V 



where 

256 




+ k- 1 + y/{K - k- 1 ) 2 + V + yj k + k- 1 - - n- 1 ) 2 + VJ . (33) 
At aurii — a22n>2 we obtain from l|32() the following approximate relation 

Pdr 1 



12 . /_3 



Pi 2aiia 22 
For n 2 S> «i and an ~ 022 we find from (|32|l that 



ni«!i. (34) 



,2 



^ w 0.8^^/^20^. (35) 

Pi a 22 

Since the scattering lengths in two-component systems are of the same order and the 
gas parameters n^a^ can be as much as 10 the quantity pdr/pi can be of order of 

10- 2 . 

Eq. H27J1 shows that the "drag density" decreases under increasing the 
temperature. Evaluating the temperature dependence of pd r numerically we conclude 
that at the temperatures of order of the interaction energy this decreasing is not 
crucial. As an example, the dependence of pdr on the temperature at n\ = = n 
and 711 = 722 = 7 is shown in Fig. Results, presented in Fig. correspond to 
?7 = 0.5, but, actually, temperature dependence of pd r is not very sensitive to this 
parameter. 

In conclusion, let us discuss to which extend the results, obtained in this section, 
can be applied to Bose gases confined in a trap. In nonuniform systems the spectrum 
of elementary excitations is modified ^3 EHl 1161 ITTj , but for healing lengths much 
smaller than the linear size of the Bose clouds the spectrum of excitations at wave 
vectors of order or higher than inverse healing length is well described by the quasi 
uniform approximation. The main contribution to the sum in Eq. (|27|l comes from 
the excitations with such wave vectors. Therefore, the drag effect in nonuniform 
system can be described by the same equations, as in the uniform case with the only 
modification that the quantities n\ and ri2, and, correspondingly, pi, p„j, pd r and ji 
in Eqs. I|25|l - (|29|l are understood as functions of coordinates. 
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Figure 1. Dependence of the "drag density" on the temperature. 



3. Model of BEC qubit with external drag force 



Let our two-component system is confined in a toroidal trap and the Bose cloud 
of component 1 (drag component) is situated inside and overlaps with the Bose 
clouds of component 2 (drive component). Such a situation can be realized if 
1 712 1 < miri (71 1,722)- Then the confining potential is deformed to cut the drag 
component in two clouds of a half-torus shape (separated by two Josephson links), but 
the toroidal shape of the Bose cloud of the drive component is survived under such a 
deformation (see Fig. |2J). In what follows we imply that circumference L = 2nR t of 
the toroidal clouds is much larger than radiuses of their cross-sections rt\ and rti- 

Rotating this trap one can excite a superflow in the drive component. After the 
rotation be switched off there will be a circulating superflow in the drive component 
and no superflow in the drag component (at negligible small Josephson coupling). The 
superfluid velocity of the drive component satisfies the Feynman-Onsager quantization 
conditions and its value is proportional to the vorticity parameter N v for the drive 
component: 

v 2 = — (36) 

m 2 R t 

(N v is integer). In we, for simplicity, neglect the effect caused by variation of rti- 
Since j\ = 0, the phase gradient V<pi should be nonzero to compensate the drag 
effect. In polar coordinates the <j> component of the phase gradient is given by the 
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Figure 2. Schematic shapes of Bose clouds for the drag (top figure) and drive 
(bottom figure) components. In system considered the drag component overlaps 
with the drive component. 



relation 

N 

(y<Pih = -IT 1 /*-, (37) 
tit 

where 

/* = ^ ^ (38) 

m 2 Pl - Pnl ~ Pdr 

is the drag factor. We imply that r t \ and r t 2 are much higher than healing lengths that 
allows to describe the drag effect in quasi-uniform approximation. For definiteness, 
we specify the case of pi <C pi and p 2 ~ const in the overlapping region. In such a 
case one can neglect space dependence the drag factor (see Eq. J2HJ). 

At nonzero Josephson coupling the current j% can be nonzero, but it cannot 
exceed maximum Josephson current J m . Relation 137fl remains approximately correct 
at nonzero Josephson coupling, if a strong inequality J rn <C Tipi/ (rnR t ) is satisfied. 
Here we specify just such a case. 

The drag force can be considered as external "magnetic field" applied to the 
drag component. The latter component in the geometry considered forms a superfluid 
circuit with two Josephson links. The situation described is analogous to the one that 
takes place in superconductive Josephson qubits ^ 111) and our system can be treated 
as a Bose qubit controlled by the drag force. 

We will formulate the model of " dragged" Bose qubit using an approach similar 
to one of Ref . [3] . 
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Let us present the Bose field operators for the component 1 in the two mode 
approximation: 

*i(r,t) = oi(*)*l(r-ri). 

(39) 



l=L,R 



*f(r,t)= ^ a+(i)*r(r-r ; ), 



l=LM. 



where o^, R -. and a L ^ are the operators of creation and annihilation of bosons in 
the condensate confined in the left(right) half-torus, and ^l, ^r are two almost 
orthogonal local mode functions 



[cPrVt(T)9i,(T)*8w, l,l' = L,R 



that describe the condensate in the left and right traps 18 . 

Substituting (|39|) into Hamiltonian (JIJ , we obtain the following expression for the 
parts of the Hamiltonian that depends on the operators Oj , a;: 



H a = ^2 (Kiafai + Xiaf a z + a ; a ; ) + (Ja^a R + J*a^a L ). 



l=L,R 



Here 



d A r^*, 



.—V 2 + V tr +712*2*2 
2m 



J 



2m 



(40) 

(41) 
(42) 
(43) 



Functions 9l and contain the phase factors e ,!P£ W and e lifR ^ with the phases 
satisfying Eq. I|37|) . Taking these factors into account, we choose the following basis 
for the one mode wave functions 



*£(«)(') = l*L(«)(r)|exp [»0i(B)(r)4*.] , 



where 



-JV»/. 



dr 



(44) 



(45) 



and 4>l,4>r are the polar angles counted from the centers of L and R Bose cloud, 
correspondingly (see Fig. |3J. The angles (f>L(R)( r ), defined as shown in Fig. satisfy 
the conditions 

4>r{?a) - 4>l{ya) = 4>l(vb) - 4>R.{rB) = 7T, (46) 

where and r# are the radius- vectors of Josephson links. 

Substituting (|44|) into Eq. 1)430. using Eq. I|46|) and taking into account, that 
the functions and *_r overlap in a small vicinity of A and B links, we obtain the 
following expression for the Josephson coupling parameter: 

J =(J A + Jb) cos(7r$ dr ) + i(J A - Jb) sin(7r$ dr ), (47) 

where 



J 



A(B) 



d A r 



V A (l 



2m 



V|* L |V|* fl |+y tr |* L ||*fl| 



(48) 




Figure 3. Left (L) and right(R) half-torus of the drag component, separated by 
Josephson links A and B. 



Here Va and V R are the areas of overlapping of two one mode functions at links A 
and B, correspondingly. 

Considering the Hilbert space in which the total number operator 

N = a^a c + a\a n (49) 

is a conservative quantity (N = N) we present the Hamiltonian (|40|l in the following 
form 

H a = E c (n RL + n g ) 2 + (Ja+a R + h.c.) + const, (50) 

where 

n RL = a+ « aR - ^ (51) 
is the number difference operator, 

E c = X R + X L (52) 
is the interaction energy, and the quantity 
1 

2E L 

describes an asymmetry of L and R half-torus. 

In what follows we imply that the system is in the Fock regime (| J\N <§; E c ) and 
use the number representation 

\ v i N N ^ 

\nRL) = \n R ,n L ) = \ — + n RL , — - n RL ). 



9 = — \Kb-K l + {N-1){\ r -\ l )} (53) 
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In this basis the first term in (|50|l is diagonal. The second term in (|50|) can be 
considered as small nondiagonal correction. But if n g is biased near one of degeneracy 
points 

!M+- for even N , , 

2 (54) 
M for odd N 

(where M is integer and \M\ < N/2) the second term in (I5U[) results in a strong mixing 
of two lowest states (| f) = Wdeg + 1/2) and | [) = \nd eg — 1/2)) and the low energy 
dynamics of the system can be described by a pseudospin Hamiltonian 

H eff = — 2° x ~ ~~2 ay ~ ~2° z ' ^ ' 

where 

** = IT>UI + IT>UI 

^ = -i(it>ai-iT>ai) (56) 
*. = it>ui-u>ui 

are the Pauli operators, and 

tt x = -{J a + Jb)tJ(N + l) 2 - n 2 deg cos(7r$ dr ), 



Sly = -{J A - J B )^j{N + l) 2 - n 2 deg sin(7r$ dr ), ( 57 ) 

Q z = 2E c (n g - rideg) 

are the components of the pseudomagnetic field. 

If the quantities n g , J a and Jb depend on time and the conditions N(Ja + Jb) *C 
E c , \n g — ridegl <C 1 are satisfied, our system can be used as qubit. In difference with 
the situation, considered in 3 , our effective qubit Hamiltonian l)55|l contains the a y 
term. The appearance of this term in the Hamiltonian is caused by the drag effect. If 
the ratio J a/ Jb depends on time and the value of 2$,jr is not integer this term cannot 
be excluded from the Hamiltonian by gauge transformation. 



4. Detection of Berry phase in the Bose qubit 

In experiments one can control the parameters n g , J a and Jb independently. It allows 
to realize an adiabatic cyclic evolution of the parameters of the Hamiltonian l|55|) , for 
which vector Q. subtends nonzero solid angle at the origin. Such an evolution results 
in appearance of a geometrical phase of the wave function known as Berry phase |19 | . 
For the system (|55H the value of this phase is the half of the solid angle subtending by 
the contour of evolution of SI. 

In this section we propose a scheme for generating Berry phase in the Bose qubit 
and show that in our scheme this phase is proportional to the value of &dr, and, 
consequently, to the drag factor. 

Let initially J a = Jb = and n g > rideg- The systems relaxes to the ground 
state | T). After that one switches on the equal Josephson coupling at links A 
and B {J a = Jb = Jin) and sets the asymmetry parameter n g equals to rideg- 
The initial state of the qubit is the superposition of two instantaneous eigenstates 
K) = (I T) + I l))/V2 and \e b ) = (| f> - I i))/V2: 

|^in) = 4=(|e„) + \e b )). (58) 
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The expectation value of the relative number operator in the initial state 

{^in\fl RL \^ m ) = 1/2. 

At the starting point the pseudomagnetic field is directed along the (-1,0,0) axes. 
Adiabatic evolution includes 4 stages. 

At the 1-st stage the Josephson coupling Jb is switched off. After this stage 
pseudomagnetic field becomes aligned along the (— cos7r<I>d r , — sin7r$(jr, 0) direction. 

At the 2-nd stage the coupling J a is switched off and simultaneously the value of 
n g is deviated from ndeg- After this stage ft is directed along the (0,0,1) axis. 

At the 3-d stage the value of n g is returned to the same degeneracy point and the 
coupling Jb is switched on. 17 becomes directed in (— cosTr^r, simt^dn 0) direction. 

At the 4-th stage the coupling J a is switched on and vector ft returns to its 
starting value. 

After such an evolution the system arrives at the state 

\^ m ) = (e i5 ^\e a ) + e^|e b » , (59) 

where 7 = ir$dr is the Berry phase equals to half of the solid angle subtending by J~2 
and d a , 5b are the dynamical phases. 

Dynamical phases depends on time of evolution, while the Berry phase does not. 
Elimination of the dynamical phases is performed by swapping the eigenstates (ir- 
transformation) and repeating the same cycle of evolution in a reverse direction |10j . 

The 7r-transformation can be done if one switches of the Josephson coupling 
and deviates the asymmetry parameter n g from rid eg value during the time interval 
t n = 2TrE c (n e ieg — n g )/h. After the 7r-transformation the state of the qubit becomes 

= --j=(e iS *+^\e b ) + e iS ^\e a }). (60) 

After the cyclic evolution in the reverse direction the state becomes 

= - A e *{*.+*)( e *7| e6 ) +e-^\e a )). (61) 

The expectation value of the difference number operator in the final state (I61|l 
differs from this value in the initial state l|58|) to the quantity that depends on the 
Berry phase and, consequently, on $>dr'- 

{ipf\n RL Hf) = ^ - sin 2 (2 7 ) = - - sin 2 (27r$ dr ). (62) 

Thus, the value of the drag can be determined from measurements of relative number 
of atoms in left and right condensates. 

5. Conclusions 

We derived a microscopic model of nondissipative drag in a two-component superfluid 
Bose gas and proposed possible implementation of a Bose qubit in which the drag 
effect is used to control quantum states. We have described how the Berry phase can 
be detected in such a qubit and found that this phase is proportional to the drag 
factor. 

Main idea of our proposal is based on a correspondence between the drag force and 
a vector potential. One can found from Eq. I|24|l that the effective vector potential, 
applied to the drag component, is equal to A e // = hfdr^^drv (in units of e = c = 1), 
where (fidrv is the phase of the drive component, and fd r is the drag factor, which 



Berry phase caused by nondissipative drag of superflow in a Bose qubit 



14 



is of order of the square root of the gas parameter. In a trap of a ring geometry 
the effective vector potential can be associated with an effective flux of external field 
& e ff — ZKhfdrNy, where N v is the vorticity of the drive component. It allows to 
realize the Bose charge qubit, analogous to superconductive qubit, considered in 
Note that while $ e ff is quantized, for fd r <C 1 it can be treated as a continuous one. 

In conclusion we would mention another possibility that opens the drag effect. 
Tuning the effective flux $ e // to the ^ff ~ n ^ va l ue an d increasing the Josephson 
coupling one can achieved a regime in which the states with different vorticities of the 
drag component are mixed. Such a system can be considered as a Bose counterpart 
of superconductive Josephson flux qubit. 
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